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Abstract — Blind Null Space Learning (BNSL) [1] has recently 
been proposed for fast and accurate learning of the null- 
space associated with the channel matrix between a secondary 
transmitter and a primary receiver. In this paper we propose 
a channel tracking enhancement of the algorithm, namely the 
Blind Null Space TY-acking (BNST) algorithm that allows trans- 
mission of information to the Secondary Receiver (SR) while 
simultaneously learning the null-space of the time-varying target 
channel. Specifically, the enhanced algorithm initially performs a 
BNSL sweep in order to acquire the null space. Then, it performs 
modified Jacobi rotations such that the induced interference to 
the primary receiver is kept lower than a given threshold Prh 
with probability p while information is transmitted to the SR 
simultaneously. We present simulation results indicating that 
the proposed approach has strictly better performance over the 
BNSL algorithm for channels with independent Rayleigh fading 
with a small Doppler frequency. 

I. Introduction 

MIMO communication has reformed the way we think 
about wireless communications, in particular, MIMO can 
be revolutionary in the paradigm of the underlay cognitive 
radios. The motivation of the underlay cognitive radios is to 
allow cognitive users to operate simultaneously in the same 
frequency band as primary users while causing them minimal 
interference ([5], [6]). 

The starting point of this work is the recent publications 
from Noam and Goldsmith ([1], [2]) where the authors provide 
a null space learning algorithm with which a MIMO secondary 
transmitter (SU-Tx) can learn the null space of the interference 
channel in order not to interfere with the primary receiver (PU- 
Rx). The algorithm has nice convergence properties, it can 
work with only energy measurements, and it is independent 
of the modulation and other transmission parameters of the 
primary system. 

Typical channel estimation problems are divided into two 
categories: channel acquisition and channel tracking ([4]). 
Channel acquisition refers to the problem of leaming the 
channel without any prior knowledge. Channel tracking refers 
to the problem of continually updating the channel estimate, 
using a previous estimate. A solution to the channel acqui- 
sition problem is provided in [3] under the assumptions of 
reciprocity of the channels, TDD transmission to both the 
primary and secondary systems, and a transmission frame 
structure of the secondary system that is synchronized with 
that of the primary system. 



In the case of time-varying channels, obtaining the null- 
space only once (i.e. null space acquisition), as it is performed 
in [1], is only the first step; a constant tracking algorithm 
that updates the null space continuously is needed. Unlike 
channel acquisition, channel tracking needs to be combined 
with simultaneous transmission of information since it is not 
a one time event, otherwise there will be a non-negligible 
reduction in data rate. Therefore, the channel tracking algo- 
rithm must ensure that the interference induced to the primary 
network remains low with high probability at all times. We 
introduce such a tracking algorithm to learn the channel and 
send information at the same time to the SR. 

The remainder of this paper is organized as follows: In 
Subsection II-A we present the channel model, in II-B we de- 
scribe briefly the BNSL algorithm and in II-C we indicate the 
performance metric that is used to test the tracking algorithm. 
Then, Section III presents simulation results that indicate that 
the null-space of a MIMO Rayleigh fading channel varies 
faster than the channel coherence time indicates. Section IV 
presents the Bhnd Null Space Tracking algorithm (BNST) 
and compares its performance with the BNSL algorithm. 
Moreover, in Section V we describe and analyze a general 
transmission scheme of information that can be used along 
with the BNST algorithm. Simulation results of the BER 
performance of this scheme are presented in Section VII. 
Finally, Section VIII concludes the paper and discusses future 
research directions. 

II. System Model 

A. Channel Model 

Consider the cognitive radio scenario depicted in Figure 1 
in which a Secondary User transmitter (SU-Tx) is allowed to 
coexist with a Primary User receiver (PU-Rx) as long as the 
interference it inflicts to him is below some threshold. Denote 
as (A^t, Ar)-system a system with Nt antennas in the SU-Tx 
and Nr antennas in the PU-Rx. To mitigate the interference 
to the PU-Rx, the SU-Tx would like to learn the interference 
channel matrix between the SU-Tx and the PU-Rx denoted 
as Hi2(t) = [hj{t)\ e C^'-^-^*, where hij{t) represents die 
channel between the SU-Tx's j-th antenna and the PU-Rx's 
«-th antenna. 

Denote as the symbol duration, as the carrier fre- 
quency, as T(t) € ([^NtxiNt-Nr) jjjg precoding matrix that the 
SU-Tx uses and as H22 {t) the channel matrix between SU-Tx 
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Fig. 1: The addressed cognitive radio scheme. 



and SU-Rx. We make the assumption of Wide Sense Stationary 
Uncorrelated Scattering (WSSUS) and that the channels hij{t) 
have independent Rayleigh fading with maximum doppler fre- 
quency Fd. This channel represents the most difficult channel 
to track since there is no correlation among the hij (t) and no 
line-of-sight (LOS) component. We also assume that between 
one cycle and another, the variation that the SU-Tx is making 
to its signal is the dominant factor that varies the PU-Rx's 
SINR. Lastly, throughout this paper we consider only flat 
fading channels. 

B. Blind Null Space Leaving Algorithm 

We now present the BNSL algorithm, for null space acqui- 
sition, from which we derive the nuU space tracking algorithm 
(BNST) tiiat is presented in Section IV. 

Assume H12 is constant and denote G = Hf2Hi2, whose 
null-space is the same as H12. The leaning process is carried 
out in cycles during which the SU-Tx transmits a learning 
signal that affects the PU-Rx interference level. We assume 
that the primary system is using a transmission scheme which 
adapts the transmitted signal according to the SINR at the 
PU-Rx. We further assume that the SU-Tx can sense whether 
the interference it inflicts to the PU-Rx increases or decreases 
from listening the PU-Tx transmitted signal or control channel. 
Denote as Tpg the amount of time needed from the time 
that the SU-Tx starts transmitting a learning message until 
it measures the effect of its learning signal on the PU-Tx's 
transmit power. Using this knowledge, the SU-Tx tries to "stir" 
the beam forming pattem of its transmissions such that it 
causes minimal interference to the PU-Rx, i.e. it transmits in 
the null space of H12. 

The learning of the null-space is carried out in sweeps. 
Each sweep consists of Nt{Nf — l)/2 learning stages where 
each one consists of K transmission cycles of length of Tp b 
msec. The SU-Tx sends constant learning signals inside each 
transmission cycle, denoted as x(n) where n denotes the 
index of the current transmission cycle. The heart of the 
"rotation" mechanism is the Jacobi Eigenvalue Decomposi- 
tion. Specifically, the Jacobi technique obtains the eigenvalue 
decomposition of the Hermitian matrix G via a series of 2- 
dimensional rotations that eliminate two off-diagonal elements 
of G in each leaming stage. The algorithm begins by setting 
Aq = G and then performs the following rotation operations 
Afe+i = B.i,,m,{0A)-^k^l,„,,{OA) where Rr„„,(^,</>) 



is an A^t X Nt rotation matrix who is equal to the identity 
matrix except for its {hJk) and {mk,mk) entries, which are 
equal to cos{9) and its {lk,Tnk) and {mk,lk) entries, which 
are equal to e^*'* sin(6') and — e""*^ sin(^) respectively. The 
rotation angles 6,(j) are chosen to zero the ilk^mk) entry 
of Afe, i.e., [Afejjj^ The integers Ik^ruk are the indices 
of the off-diagonal entries that are eliminated in the fc-th 
step and are chosen such that the off-diagonal elements are 
eUminated periodically, where each entry is ehminated once in 
each leaming stage [ 1 ] . Noam and Goldsmith proposed a blind 
technique (without observing the matrices A^^k = 0, 1, ... ) 
in which the SU-Tx obtains 6 and with only two binary line 
searches. In other words, the SU-Tx in each learning stage k 
calculates a pre-coding matrix Wfc such that W^GWfe = A^. 
To do so, it samples the ^ — </> plane and gets feedback from 
sensing the transmit power of the PU-Tx. After Nt{Nt - l)/2 
leaming stages , which constitute one Jacobi sweep, the SU-Tx 
has the first estimate of the eigenspace. 

C. Performance Metric 

In order to compare the different null space leaming meth- 
ods we need to introduce a natural performance metric. The 
metric that the primary system is concemed about is the 
maximum interference caused by the SU-Tx. Denote as x(t) 
the message that the latter wants to send. Then, the normalized 
interference is upper bounded by the spectral norm of the 
matrix Hi2(i)T(t) as: 

||Hi2(/.)Ti/)x(t)||2 



|Hi2WT(t)|| 



max ,, , 

||x(t)||2^0 ||x(i)||2 
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We say that the null-space tracking is successful if 
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that is, if p percent of the time the interference achieved with 
the use of the precoding matrix T{t) is at least Pth dB less 
than if we use as T{t) = iNtxNf 

III. Null space coherence time 

In this section we motivate the need for using the BNST al- 
gorithm based on the observation that even in slowly Rayleigh 
fading MIMO channels, the null space changes much faster 
than the channel coherence time associated with each of the 
independent channels hij{t) indicates. 

The coherence time T,. is a widely accepted characteristic 
of a channel that is used to quantify how fast it changes. We 
assume that the correlations between channel coefficients at 
different times depends only on the time difference At. We 
define the normahzed autocorrelation of a SISO channel h{t) 
as (chapter 3.3 [4]) 

E{h{t)h{t + At)*} 



PAt 



Emm 



(2) 



Then the X% coherence time, denoted as tF>, is defined 
as the value At such that p^t = f^. i.e. Pj,[x] = For 

example, in Clarke's model [4], j'J^'^l = jq^^. 



To study the relation between the Hi2(t)'s null space time 
variations and the coherence time of its entries, we turn to 
simulations. Consider a (3, l)-system where each hij{t) is 
an independent Rayleigh fading channels (i.e. the entries of 
Hi2(t) are independent of each other) with Doppler Frequency 
of Fj, — 6.48 Hz (e.g. relative speed of 10 Km/h and 
carrier frequency of Fc — 700 Mhz) using the Clarke's 
model with 40 multipath components. Figure 2 depicts the 
amplitude of the normalized autocorrelation IpAtI and the 
relative interference reduction dMii^t) i.e. 



|Hi2(t + At)AA(Hi2(0)|| 
l|Hi2(t + At)|| 



dMi{At) = lOlogio E 



(3) 



The latter calculates the average minimum decrease of in- 
terference inflicted to the PU-Rx if the current channel is 
Hi2(i + At) whereas the SU-Tx uses the null space of 
an outdated channel, i.e. JV(H.i2{t)). Note that dMii^t) 
increases surprisingly fast. For example, dMi > — 15dB for 
pAt < 0.95! This means that even if the SU-Tx obtains 
the null space of the channel instantly at time t, it will 
have to update it after T^^-^^^ seconds, if it needs to cause 
a minimum of —15 dB interference to the primary system. 
The intuition behind this observation is that Hi2(i) is a time- 
varying matrix whose elements change approximately every 
Tc seconds. The null space is a quantity that depends on all 
the elements of the matrix and therefore it is likely to change 
much faster. This observation was verified using a different 
Rayleigh fading generator ([7]) and with different numbers of 
multipath components (from 5 up to 100). 
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Fig. 2: (3, l)-system. The left vertical axes represents 
dMii^t) and the right vertical axes pAt- Horizontal axes 
represents the normahzed time as a a fraction of 



IV. Blind Null Space Tracking 

In this section, we summarize the Blind Null Space Tracking 
(BNST) algorithm and present simulation results which show 
that in the case of slow-varying MIMO channels it successfully 
tracks the channel. We postpone the discussion on how the 
SU-Tx performs simultaneous transmissions while it tracks 
the null space in Section V. 



Denote as T the time needed for a Jacobi sweep to finish. 
Assume at time t we have an approximation W(i) of the 
eigenspace W(i) of the channel matrix. 

Consider the Modified Jacobi Sweep that is described in 
Algorithm 1 . Function NextElement returns the values Ik , nik 
for the fc-th learning stage. The input is the current estimated 
eigenspace W(t) along with the parameters 9,9max,V that 
define the search space of the line searches of the sweep. 
Specifically, Omax denotes that the SU-Tx should perform a 
line search in the interval [—dmax,Smax]- S is the value of 6 
used in the LineSearch in Line 4 of Algorithm 1 and rj is the 
step size of the binary search which may vary as a function 

of Qraax- 

The channel tracking works as follows. Firstly, the SU-Tx 
performs one complete sweep of the BNSL algorithm, i.e. it 
uses W(i) = liqt>^Nt,0 = ^,Omax = f , = 0.1. Thcsc are 
the parameters used in [1]. This sweep creates the first estimate 
W(< + T) of the eigenspace and completes the channel 
acquisition phase. Then, the SU-Tx transmits information 
using as the precoding matrix T = [wi, wi, . . . , watj-at^], 
where Wj are the columns of the matrix W(f + T), i.e. the 
SU pre-multiplies its transmitted signal t + T by the matrix 
T. The modified BNSL sweep always returns the W(t + T) 
such that the singular values that correspond to its column 
are in a decreasing order, i.e. let a-i be the singular value that 
corresponds to w;, then <7i < Ui+i, 1 < i < Nf. 

The SU-Tx continuously senses the transmit power of the 
PU-Rx. When it senses that the transmit power has increased 
more than a predefined threshold, then it updates the null 
space's estimate. In our current simulations we make the as- 
sumption that the SU-Tx knows ' the quantity ||Hi2(t)T||rfB 
and it performs an adaptation whenever | |Hi2(t)T| j^is < Pxr 
where Pxr = —20. The intuition behind this rule is that 
the SU-Tx should not wait until the current eigenspace has 
changed significantly since in that case the matrix W(t + T) 
cannot be of much help in the subsequent sweep. 

In order to update fast, accurately and with minimal interfer- 
ence caused to the primary system, the algorithm performs a 
local search over the space of the angles. This local search 
should be performed for the following two reasons. Firstly, 
each line search can be much faster, assuming that the SU-Tx 
keeps the same parameter rj. Secondly, given that the SU- 
Tx has a good estimate of the null space, there exists an 
explicit trade-off between how large O^ax should be and the 
inflicted interference to the primary system . To be precise, it 
is known that if the estimate is very close to the real null space 
then 9^0 ([1]), i.e. there is no need of a "large rotation". 
Therefore, the SU-Tx should not search for the optimal 6 over 
a large interval around since it knows that it has a good 
approximation already. A similar reasoning explains why 9 
is an important parameter of the BNST; the SU-Tx should 
perform the first line search (Line 4 in Algorithm 1) while it 
is transmitting close to the approximated null space. If the SU- 



' In [1], it is shown how can the SU autonomously estimate the interference 
reduction without cooperation from the primary user 



Algorithm 1: Modified Jacobi Sweep 
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Output: W(t + T) 
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Tx has a good estimate of the eigenspace, a small means that 
the line search in line 4 of the modified rotation algorithm is 
performed around the null space of the channel. Then, the SU- 
Tx manages to adapt and still keep the inflicted interference 
in low levels. 

Note that if the estimate is severely outdated, either because 
the SU-Tx waited too long to perform another sweep, or 
because the doppler frequency is large, then the SU-Tx should 
search for the best 9 over the whole interval [— 7r/2, 7r/2]. In 
our simulations O^ax was and ^ for small ([1 — 2 Hz]) and 
large (> 2 Hz) doppler frequencies respectively. A complete 
sweep uses Omax = f- The parameters 6, 9 max should be 
chosen as a function of the doppler frequency Fd- These can 
be predefined quantities that the SU-Tx chooses after it has 
estimated the doppler frequency of the channel. 

In Figure 3 we show an example {Nt — 2, Nj. — 1, — 
66.7 /isec, Tpb — 1 msec ^) of the minimum decrease in 
interference (in dB) as a function of time when Hi2(t) is a 
time varying channel with Fd ~ 1.3 Hz. 

In Figure 4 we compare the performance of the BNSL 
algorithm against the performance of the BNST algorithm over 
a (2, l)-system. We average over 100 different channels and 
10* slots for each channel for different doppler frequencies 
Fd- The performance metrics that we use are P95, P90, Pss 
and Average Decrease of Interference, where Px represents 
the decrease of interference (in dB) in X% percentage of the 
slots compared to the system that does not use any precoding 
matrix. For instance, if Pgg — — 10 dB it means that in 90% of 
the slots the SU-Tx caused at least —10 dB less interference 
by using the proposed algorithm. Note that whenever the SU- 
Tx adapts we calculate the actual interference caused to the 
PU-Rx since the SU-Tx sends known signals defined from 
the algorithms. On the other hand, when the SU-Tx is not 
adapting we use the maximum interference as defined in (1), 
i.e. assuming the worst case scenario. 

We observe that the BNST algorithm performs significantly 
better as expected since the BNSL algorithm can be considered 

-These parameters are in agreement with the LTE current specifications. 
Any decrease in Tpg will lead to a more accurate and fast tracking. 



as a special case of the BNST algorithm. However, we 
observe that for high doppler frequencies, the BNST is not 
able to keep Pgs in low values, which means that the null- 
space changes faster than the adaptations are performed and 
therefore, even when SU-Tx performs a modified adaptation it 
causes significant interference. On the whole, we observe that 
for low doppler frequencies, even in this worst case scenario of 
independent Rayleigh Fading MIMO channels, the algorithm 
manages to keep the interference at low levels. 
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Fig. 3: Example for a (2, l)-system with Fd = 1.3 Hz. 




Fig. 4: Pg5, P90, Pss and Average decrease in interference 
(Label: Av) for the BNSL and BNST algorithms for different 
Fd- Superscript L and T mean BNSL and BNST respectively. 



V. Simultaneous Data Transmission And Tracking 

As we have already observed, even if the SU-Tx manages 
to track the channel, it will be able to transmit information 
to SU-Rx for a small number of time-slots and then it will 
need to adapt again. Therefore, the secondary system uses 
a significant amount of its transmission time only to adapt 
to the changing environment. Motivated by this observation, 
the SU-Tx uses the previously presented ideas in order to 
satisfy the interference constraints imposed by the primary 
network and then we propose a transmission scheme that 
can be employed to transmit information superimposing the 



information onto the learning signals without introducing any 
error in the tracking. 

Denote as iV = [^^^] the number of consecutive slots 
that the SU-Tx transmits the same learning vector, as r2(<) 
the information signal that SU-Tx wants to send to SU-Rx, 
and as H22(i) the channel between the SU-Tx and SU-Rx. 
During the TpB period the channels Hi2(f) and H22(t) will 
be approximately constant, i.e. 

Hi2(i) « Hi2 and H22(t) « H22 

because TpB « Tc, otherwise the tracking would be impos- 
sible as we have already discussed in Section III. 

The idea is to find a way to superimpose r2{t) so that the 
functionality of the BNST algorithm remains unaffected and 
the SU-Rx receives information. We start from the signal that 
the PU-Rx and the SU-Rx receive respectively: 

yi(i) = Hi2T(ri+r2(i)) + ni(i) (4) 
y2(i) = H22T(ri+r2(i)) + n2(t) (5) 

Denote also as yi{t) and y§(f) the received signals at the 
PU-Rx and the SU-Rx respectively when no information is 
superimposed to the learning signal, i.e. 

y? (t) = Hi2Tri + m (i) , y° {t) = HasTn + na (i) (6) 

Define as 

At/1 = Q{yi) - Q(y?) (7) 

the difference in the measurement at the PU-Rx as a result of 
the introduction of the r2{t) signal, where Q(-) is a measure- 
ment function Q : R used by the PU-Rx in order to 
estimate the interference. The SU-Tx needs to superimpose 
the information over the learning signal in a way that the 
output of the <5( ) function remains unaffected. Therefore, the 
transmission scheme depends on the fimction Q{-). We are 
going to consider two different energy measurements that can 
be used in practice. 

A. First Model of Energy Measurement 

Assume that the PU-Rx measures the interference with the 
following way: 

1 ^ 

i.e. for each time-slot it measures the power of the interfering 
received signal and then it calculates the time-average over iV 
slots. 

We make no assumption that the SU-Rx knows the current 
precoding matrix of SU-Tx nor the channel H22 between 
them. From equations (4),(6),(7) it follows that 

A2/i = Qi(yi)-0i(y?) 

1 ^ 

= i\rE(l|Hi2T(ri+r2(t)) + ni(t)||2 
-||Hi2Tri+ni(i)||2) 



A sufficient condition for Y2 (t) to not affect the Jacobi method 
is: Ayi = 0. Assume that 

r2(t)=c(f)ri (9) 

where c{t) G C is drawn according to a known discrete 
probability distribution on the support {ci, C2, • • • , cm}, and 
M is the number of different messages that can be sent. A 
reasonable assumption would be Pr[c{t) = Cj] = Vi € 
{1, 2, • • • M}, but we do not restrict this scheme only to this 
case. This assumption is the heart of the transmission scheme 
and it means that the SU-Tx uses the learning signal ri as 
the "carrier" signal of transmitting information. Consider the 
following constraints that the {c;} should satisfy: 

• E{|1 -|- c(t)p} = 1. This is sufficient in order to ensure 
that the Ayi is close to 0. Note that this can always be 
ensured if |1 + c^p = l,Vi e {1, 2, • • • M}, but still this 
is not the only way. 

• E{c(i)} 7^ —1. The necessity of this assumption will 
become evident later on. 

It follows that 

= ^ E (IIH12T [ri + r,{t)] + ni(i)||2 

-||Hi2Tri+ni(i)in 
"4-' E{||Hi2T(ri +c(i)ri)|p} - HHiaTVilp (10) 
= ||Hi2Tri||2E{|l + c(0|'}-||Hi2TVi||2 ^=^0 

where 

• w.p.l means with probability 1, 

• equation (10) follows from the strong law of large num- 
bers (LLN). 

. (6) follows from the fact that E{|1 + c(i)|2} = 1. 
Thus, Ayi w for sufficiently large N. 

Now, we need the SU-Rx to be able to decode the received 
message. In the end of the N slots, the SU-Rx can estimate 
the average of the received signal: 

1 ^ 

y2 = ]vEy2(t) (11) 
t=i 

1 ^ 

= ]7E(H22T(ri+r2(0) + n2(t)) (12) 

i=l 

"'4' H22Tri + H22Tri (E{c(<)}) = CH22Tri (13) 

where 

• equation (13) follows from applying the LLN to both 

ni(i) and r2(t) 

• C = 1+E{c(f)} is a predefined parameter known to both 
the SU-Tx and SU-Rx. 

Denote as e(A^) = CH22Tri — y2 the error introduced 
because N is finite. Then, after N slots the SU-Rx can decide 
which Ci was sent by subtracting the value ^y2 from y2(i): 

Ay2(t) = y2W-^^y2 (14) 
= (c(f) - Ci)H22'IVi + f{e{N), n2{t)) (15) 



where /(e(7V), n2(i)) = (1 + c^)^ + n2(t). Note that 
/(e(iV),n2(t)) -J> n2(t) for -J> oo and that T is not in 
the null space of H22 with probabihty 1([2]). 

Thus, the SU-Rx can decide amongst which Ci was sent 
in each time slot. More specifically, consider the set of 
hypothesis {Hi} such that Hi = |4y2- Since Ci have equal 
prior probability and assuming n2 (t) is white Gaussian noise, 
the probability of error will be minimized if we foUow the 
following decision rule: 

data(t) = argmini<i<M{||Ay2(t) - Hi\\} (16) 

The main assumption that we have made in the above 
discussion is that N is large enough in order the LLN to 
hold. Even if this is not true, we can ensure our performance 
requirement is met using Enumerative Coding (Section VI). 

B. Example for a binary alphabet (M = 2) 

Assume that the SU-Tx is using a 2-constellation signal on 
top of ri 



In this case, the proposed scheme can be modified appropri- 
ately. Again, we do not assume that the SU-Rx knows the 
current precoding matrix T nor the channel H22. Therefore, 

Ayi = g2(yi)-Q2(y?) 



r2(t) = 



ciTi, SU-Tx transmits 1 with prob. 0.5 
C2ri, SU-Tx transmits 2 with prob. 0.5 



where ci and C2 are complex numbers. The constraints that 
the ci and C2 should satisfy are: 

. |l + Ci|2 + |l + C2p = 2. 

• ci+C2 7^ -2. 
If we restrict to the case that ci = e-'^^ and C2 = e~^^° then 
both constraints are satisfied for 9o = Notice that in this 
specific case we always get: 

||Hi2Tri(l+e^'^'')||2 = ||Hi2'IVi(l+e-^''°)||' = HHiaTnH^ 
Then, C = 1 + cos(0o) and 

_ f T+f^sIeo)^2(*) + n2(t), if 1 was transmitted 
I Tq^^^y2(t) + n2{t), if 2 was transmitted 

The two assumptions Hi and H2 of the decoder are: 



Hi 



1 + cos(6'o) 
The decision rule is 



Mt), H2 



1 + cos(6'o) 



y2{t) (17) 



data(t) 



1, if\\Ay2{t)-Hi\\<\\Ay2{t)-H2\ 

2, other 



C. Second Model of Energy Measurement 

It is possible that the PU-Rx uses the more reUable energy 
measurements of the form: 



1 ^ 

Q2(yi) = ii^EyiWii 



(18) 



,^|^(Hi2T(ri+r2(t)) + ni(i))|p 
t=i 

-|li\rE(Hi2TVi+ni(t))|p 



(19) 



Assume again as before that T2{t) = c{t)ri and consider the 
following constraints on c^: 

. |l+E{c(i)}P = l 

• Note that the constraint E{c(t)} 7^ —1 is inomediately 
satisfied. 

Starting from (19) and using the LLN as before: 

w.p.l 



|Hi2TVi(l +E{c(0})|r - IIHialVill^ = 



The SU-Rx in this case works exactly the same as before so 
we do not repeat the procedure. Note, that the only difference 
is that the set {cj} needs to satisfy different constraints. 

VI. Enumerative Coding 

In the above scheme we made the assumption that is 
large enough for the LLN to hold. To be precise, if we can 
ensure that 



1 ^ 



(20) 



for all possible transmitted sequences c(l), c(2), • • • ,c{N), 
where 5 is a predefined constant known to both the SU-Tx 
and SU-Rx, then the LLN assumption is not needed at all. 
Assume the secondary system fixes S, then it needs a coding 
strategy that maps an iV-tuple of symbols to another A/'-tuple 
that satisfies the constraint. 

A simple and efficient method that can ensure this constraint 
has been treated by Cover ([8]). Specifically, without loss of 
generality, assume that the SU-Tx wants to send symbols from 
the alphabet X = {1,2, ...,M} according to the uniform 
distribution. This assumption models well the distribution of 
the actual data that are sent in a wireless network. We need 
some definitions from the information theoretic literature. 
Define the type Px of a sequence x = x" as the relative 
proportion of occurrences of each symbol. Also, define as 
the type class of the probability distribution P, denoted as 
r(P) the set of sequences of length N that have type Px, i.e. 

r(P) = {x e A"^ : Px = P}. 

For simplicity we assume that is an integer. 

Equation (20) is satisfied if the transmitted sequence 
c(l),c(2),...,c(7V) has P^ = ijj^jj,--- , m), or in other 
words, if the secondary system uses only those sequences that 
belong in the type class T(Pc). We know that 



\T{Pc 



N 

N_ N_ N_ 
M ' M ^ ■ ■ ■ ' M 



(21) 



Therefore, the SU-Tx can send only |T(Pc)| different A^- 
sequences from the available sequences. For example, 
for M = 2 and iV = 16 we get that [log2(|r(Pc)|)J = 13. 
For M=2 and N=32, [log2(|T(Pc)|)J = 29 and for M=4 and 
N=16, Llog2(|T(P,)|)J =25. 

The encoding process is the following: Assume that the 
SU-Tx wants to transmit the sequence a;", where n — 
[log2(|T(Pc)|)J . Then, it treats the sequence x" as the index 
of an iV-sequence that belongs to the type class r(Pc). 
Thus, it maps the index to a sequence and then transmits 
the sequence. The SU-Rx performs the reverse procedure. 
Both encoding and decoding procedures are more thoroughly 
explained in [8] and can be easily performed with low com- 
plexity. 

VII. Simulation Results 

In this section we test the effectiveness of the proposed 
scheme in its most generality, i.e. outside the BNST algorithm. 
Specifically, we generate 10"^ random matrices H12, H22 and 
a random unitary matrix T. Then, for each simulation scenario 
we generate 10^ bits that are superimposed on a vector ri that 
was created by a random Jacobi rotation. We use = 16, 
M = 2 (assuming TpB = 1 msec and = 66.7/isec) and 
6*0 = ^ (See Section V-B). Pg is the average probability of 
symbol error over all scenarios. In Figure 5 we plot the P, 
as a function of the SNR at the the SU-Rx using the scheme 
presented in section V-B. 
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Fig. 5: Ps as a function of SNR at the SU-Rx. 

The average value of Aj/i is 0.13 dB, which is an insignifi- 
cant increase of the interference. This means that the learning 
procedure remains unaffected by this transmission scheme. 
From Figure 5 we observe that in the case of an (3, l)-system 
and (2, l)-system and M = 2, for SNR> 3.5 dB or SNR> 4.5 
dB respectively, the probability of error is less than 10^"^. We 
observe that more antennas in the SU-Tx leads to a more 
robust system. Similar results are shown also for M = 4. 

VIII. Conclusions - Future Work 

The starting point of this work was the observation that 
the null space of a time-varying channel changes significantly 
faster than the coherence time of the channel. This motivated 
the need for a null space tracking algorithm, based on the ideas 
of the BNSL algorithm, which learns the null space without 



inflicting excessive interference and transmits information 
simultaneously over the secondary system. We proposed the 
BNST algorithm in order to track the null space variations 
and showed with simulations that the proposed algorithm can 
actually enhance significantly the performance of the BNSL 
algorithm for low doppler frequencies. Lastly, we presented 
a transmission scheme that coexists with the learning process 
by superimposing the information on the learning signal and 
demonstrated its validity through simulations. 

No matter which is the algorithm for null-space acquisition 
and/or tracking we observed through simulations that the null 
space under independent Rayleigh fading changes significantly 
faster than the coherence time. More research is needed in 
order to understand how fast the null space varies. Also, the 
Rayleigh fading MIMO channel with independent coefficients 
is obviously the worst case scenario since there is no LOS 
component. We could test the results of proposed approach in 
difference channel scenarios, for example when the channel 
has Rician fading. 

Last but not least, in this work we used the vector ri, i.e. 
the current learning signal, as the carrier of the information we 
need to send. This gives us only one spatial degree of freedom. 
Is it possible to propose a similar algorithm with which we are 
transmitting using two or more spatial degrees of freedom and 
still we do not affect the learning process? These directions 
are going to be under investigation in our future work. 
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